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ABSTRACT 

The idea that non steady accretion during the embedded phase of protostar evolution can 
produce the observed luminosity spread in the Herzsprung-Russell diagram (HRD) of young 
clusters has recently been called into question. Observations of Fu Ori, for instance, suggest 
an expansion of the star during strong accretion events whereas the luminosity spread implies a 
contraction of the accreting objects, decreasing their radiating surface. In this paper, we present 
a global scenario based on calculations coupling episodic accretion histories derived from numer- 
ical simulations of collapsing cloud prestellar cores of various masses and subsequent protostar 
evolution. 

Our calculations show that, assuming an initial protostar mass M; ~ lMj up , typical of the 
second Larson's core, both the luminosity spread in the HRD and the inferred properties of 
Fu Ori events (mass, radius, accretion rate) can be explained by this scenario, providing two 
conditions. First, there must be some variation within the fraction of accretion energy absorbed 
by the protostar during the accretion process. Second the range of this variation should increase 
with increasing accretion burst intensity, and thus with the initial core mass and final star mass. 
The numerical hydrodynamics simulations of collapsing cloud prestellar cores indeed show that 
the intensity of the accretion bursts correlates with the mass and initial angular momentum of 
the prestellar core. Massive prestellar cores with high initial angular momentum are found to 
produce intense bursts characteristic of Fu Ori like events. Our results thus suggest a link between 
the burst intensities and the fraction of accretion energy absorbed by the protostar, with some 
threshold in the accretion rate, of the order of lO _5 M0yr _1 , delimitating the transition from 
"cold" to "hot" accretion. Such a transition might reflect a change in the accretion geometry 
with increasing accretion rate, i.e. a transition from magnetosphcric or thin disk to thick disk 
accretion, or in the magnetospheric interaction between the star and the disk. Conversely, the 
luminosity spread can also be explained by a variation of the initial protostar mass within the 
<~ 1-5 M Jup range, although it is unclear for now whether such a spread among second Larson's 
core can be produced during the prestellar core second collapse. This unified picture confirms 
the idea that early accretion during protostar and proto-brown dwarf formation/evolution can 
explain the observed luminosity spread in young clusters without invoking any significant age 
spread, and that the concept of a well defined birthline does not apply for low-mass objects. 
Finally, we examine the impact of accretion on the determination of the IMF in young clusters. 

Subject headings: stars: formation — stars: low-mass, brown dwarfs — accretion, accretion disks 
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during the embedded phases of protostellar evolu- 
tion could explain the luminosity spread observed 
in young, a few Myr old, clusters and star forming 
regions. The same scenario also provides a natu- 
ral explanation for the unexpected lithium deple- 
tion observed in some young objects (Baraffe & 
Chabrier 2010, BC10). Given the increasing evi- 
dence that low-mass protostars accrete a substan- 
tial part of their material through short intense 
bursts of accretion (Kenyon et al. 1990; Enoch 
ct al. 2009; Evans et al. 2009), due to gravita- 
tional (Vorobyov & Basu 2005, 2006; Machida et 
al. 2010) or a combination of gravitational and 
magneto-rotational (Zhu et al. 2009) instabili- 
ties in the accretion disk, the scenario proposed 
by BCG09 is based on the process of so-called 
" episodic accretion" . Assuming that the accreting 
material radiates away essentially all its internal 
energy at the shock front, with no contribution to 
the internal heat content of the accreting object 
(the so-called "cold accretion" process, see §3), 
and taking for the initial protostar a mass range 
from 1 Mj up to 0.1 M Q , BCG09 were able to pro- 
duce a spread in luminosity in the Herzsprung- 
Russell diagram (HRD) for the final objects, at 
the end of the accretion process, i.e. after a few 
Myr, equivalent to a typical <~ 10 Myr age spread 
for non-accreting objects. Their major conclusion 
was thus that the observed HRD luminosity spread 
does not stem from an age spread, but rather from 
the impact of early accretion history on the struc- 
ture and the cooling of the protostar. This conclu- 
sion is supported by the recent observational anal- 
ysis of stars with and without disks in the Orion 
Nebula cluster by Jeffries et al. (2011). Finding 
no significant difference in the mean ages or age 
distributions inferred for stars with and without 
disks, respectively, these authors conclude that the 
observed luminosity dispersion cannot be due to a 
large age spread, but rather stems from a combi- 
nation of observational uncertainties and physical 
mechanisms. 

Conducting a study similar to BCG09, Hosokawa 
et al. (2011, HOK11) recently reexamined the im- 
pact of early accretion on the evolution of low mass 
stars (LMS) and brown dwarfs (BD). Although 
these authors find results similar to BCG09 con- 
cerning the effects of accretion on the structure 
of these objects, confirming BCG09's study, they 
reach different conclusions and claim that accre- 



tion cannot produce the observed spread in the 
coolest part of the HRD, i.e., for objects with 
T eff < 3500 K. Further contributing to the de- 
bate, Hartmann ct al. (2011, HZC11) recently 
reanalysed the spectral energy distribution of FU 
Orionis. This object undergoes short accretion 
bursts from a disk and provides an excellent labo- 
ratory to test the effects of strong accretion bursts 
on the structure of accreting objects. HZC11 high- 
light the absence of a hot boundary layer emission, 
which is best explained by significant heating of 
the protostar upper layers by the disk infalling ma- 
terial. Furthermore, HZC11 infer a rather large 
radius (<~ 5 R ) for the protostar whereas, in 
case of aforementioned cold accretion, BCG09 pre- 
dict a much more compact structure, because of 
the accretion-induced contraction of the object. 
HZC11 thus conclude that for accretion histories 
with burst rates as high as about 10 _4 MQyr _1 , as 
observed in Fu Ori objects, some fraction of the 
accreting energy must be absorbed by the proto- 
star (departing from the genuine "cold accretion" 
regime), yielding an expansion of the object, as in- 
deed obtained by BCG09's in the "hot-accretion" 
case. 

In order to clarify these controversies, we 
present in this paper a global scenario, based 
on the coupling between episodic accretion cal- 
culations obtained from collapsing cloud prestel- 
lar core simulations, and accreting protostar/BD 
evolutionary models. We show that such a sce- 
nario explains both the HRD luminosity spread 
and the expansion of Fu Ori stars, reinforcing the 
conclusions reached in BCG09. We also explain 
why HOK11 reached different conclusions than 
BCG09, even though conducting similar calcula- 
tions, which should help resolving the controversy 
raised about this issue. 

2. Accretion histories from prestellar core 
collapse calculations 

In the present calculations, the accretion rates 
on the protostar, M(t), are obtained from hydro- 
dynamics simulations of prestellar core collapse 
similar to the ones conducted in Vorobyov & Basu 
(2010, VB10) and Vorobyov (2010, V10). All the 
details of the model can be found in VB10. The 
simulations start from collapsing prestellar cores 
of mass M c with Bonnor-Ebert like density pro- 
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files, and follow the formation and evolution of 
the protostar (mimicked as a point mass gravitat- 
ing object at the centre of a sink cell of radius 
r sc = 6 AlQ ), the disk and the envelope. The 
mass accretion rate is calculated as the mass pass- 
ing through the sink cell per time step of numeri- 
cal integration (which in physical units is usually 
equal to 10-20 days and the total integration time 
can amount to 2.0 Myr). The protostar acquires 
a fraction £, ranging between 50%-90%, of the ini- 
tial prestellar core's mass, implying that a fraction 
1 — £ is assumed to have been carried away by out- 
flows)!] 

In the present calculations, the structure and 
evolution of the accreting protostar is consistently 
calculated from accreting evolution models (see 
§3). We stress that these calculations take into 
account the thermodynamic properties of the gas, 
i.e., protostellar and background irradiation, ra- 
diative cooling and viscous and shock heating (see 
VB10). We have conducted simulations, and thus 
obtained accretion histories, for initial prestel- 
lar core masses M c ranging from 0.06 Mq to 
1.85 Mq. The strongly variable accretion rates 
obtained in this type of simulations have been 
shown to reproduce the observed steep M-M* re- 
lation in class II objects (Vorobyov & Basu 2009). 
They also reproduce the observed spread for pro- 
tostar signatures both in the L^ \-T^ \ and L^ Q \- 
M c diagrams, including the very low luminosity 
objects (VELLO's) (Dunham & Vorobyov 2011). 
Even more interestingly, they provide an excellent 
match to these protostellar observed luminosity 
distributions, strongly supporting episodic accre- 
tion as the most plausible explanation for the long- 
standing luminosity problem. 

VB10 and V10 have explored the dependence 
of the results upon the initial core mass, M c , 
and rotational to gravitational energy ratio, /3 = 
E rot / E grav , within a range of values (j3 = 10~ 3 — 
3x 10~ 2 ) typical of observations of dense molecular 
clouds (Caselli et al. 2002). They show that high 



1 Our test runs have shown that decreasing r sc by a factor of 
three does not qualitatively change the accretion pattern 
but entails almost a factor 10 increase in computational 
time, making simulations of a large number of models with 
small sink cells prohibitively expensive. 

2 As mentioned in VB10, the smaller the fraction of core 
mass ending up as the protostar mass, the larger the disk 
to star ratio, favoring more disk instability and thus burst 
accretion. 



M c and high j3 favor disk fragmentation, because 
of the increase in the disk mass due essentially to 
the larger initial core mass or the larger centrifugal 
radius for the more massive or more rapidly rotat- 
ing cores. Given that the accretion burst intensity 
correlates with the disk's propensity to fragment 
(the bursts are caused by massive fragments merg- 
ing with the protostar due to gravitational torque 
interactions with the spiral arms), the accretion 
burst intensity is thus found to increase with M c 
(see e.g. Fig. 8 of VB10) and with f3 (see e.g. Figs. 
5 and 7 of VB10 and Fig. 6 of V10). It is worth 
stressing that in these simulations, only a combi- 
nation of a high prestellar core mass, M c > 0.7 M Q , 
and a high rotation rate, /3 > 5 x 10 -3 , can produce 
burst intensities > 10~ M©yr - , as inferred for Fu 
Ori type objects (HZC11, sec §1). Similar results 
have been found in 2D and 3D simulations explor- 
ing the same issue (Machida et al. 2010, Zhu et 
al. 2011). 

These prestellar core collapse simulations have 
two major limitations. First, they are limited to 
2D, though with an approximate reconstruction 
of the disk vertical structure. The main reason 
for such a limitation is that 2D simulations allow 
a numerical resolution which is presently out of 
reach with 3D simulations, an important issue to 
accurately describe the disk structure and evolu- 
tion, in particular when taking into account the 
thermodynamics of the gas. The second reason is 
that 2D simulations enable us to cover many disk 
orbital periods (with absolute integration times of 
order 1.0-2.0 Myr), which is not possible so far in 
3D, again a major issue when exploring the disk's 
long-term evolution and stability and the fate of 
the massive fragments. A first set of 3D simu- 
lations exploring the stability of the disk during 
protostar formation has recently been obtained by 
Machida et al. (2010, 2011). These simulations 
qualitatively and even semi-quantitatively agree 
with the 2D ones of Vorobyov & Basu, exhibit- 
ing the same episodic accretion phenomena, with 
comparable burst intensities. The second limita- 
tion of the present simulations is the absence of 
magnetic field. It is well known that magnetic field 
plays a crucial role in the prestellar core collapse 
and can significantly hamper disk fragmentation, 
at least in ideal MHD calculations (Banerjee & 
Pudritz 2006, Price & Bate 2007, Hennebelle & 
Teyssier 2008, Machida et al. 2008, Hennebelle et 
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al. 2011). The ideal MHD simulations of Vorobyov 
& Basu (2006) have shown that magnetic fields can 
significantly reduce the burst intensity, although 
this effect can be counterbalanced by increasing 
the initial mass and angular momentum of the 
cloud prestellar core. The importance of magnetic 
field upon the formation of the proto circumstel- 
lar disk has received support from confrontation 
with observations (Maury et al. 2010). Recent 
resistive MHD simulations, however, suggest that, 
although considerably smaller than that in an un- 
magnetized cloud, a circumstellar disk can still 
form in a strongly magnetized cloud (Inutsuka et 
al. 2010, Machida et al. 2011), even though the 
issue does not appear to be completely settled yet 
(Li et al. 2011). The orientation of the field is 
also found to have an impact, with departure of 
the field from alignment with the rotation axis 
favouring the formation of a disk (Hennebelle & 
Ciardi 2009). Therefore, although the exact burst 
characteristics obtained in the present 2D, hydro- 
dynamics calculations may be questionable from 
the very quantitative point of view, the general 
properties of the accretion histories seem to be 
reasonably sound. 

Finally, a limitation of the present calculations 
is the lack of self-consistent treatment of the radia- 
tive feedback from the accreting protostar. Our 
numerical hydrodynamics code (VB10) uses ac- 
cretion and photospheric luminosities calculated 
from the non-accreting stellar evolution models of 
D'Antona & Mazitelli (1997), which differ from 
those calculated from accreting stellar evolution 
models. Fully self-consistent simulations require 
the coupling of our hydro and stellar evolution 
codes in one real-time numerical simulation of the 
stellar and disk evolution during the main accre- 
tion phase. Work to produce such self-consistent 
models is under progress (Vorobyov, Baraffe & 
Chabrier, in prep.). 

3. Accreting models for the protostar and 
brown dwarf evolution 

The general energy equation for an accreting 
object reads (Chabrier et al. 2007a, BCG09): 

, GMM GMM 

Ltot = (l-a)e — hae — - — 

K K 

- J M T{ ^- m ^ )t}dm + LD 



(1) 

The first term of the r.h.s. of cqn.(l) corresponds 
to the radiated accretion luminosity, L acc , the sec- 
ond term to the fraction of accretion kinetic energy 
absorbed by the protostar; the third term corre- 
sponds to the protostar intrinsic entropy rate at 
constant mass shell q = m/M plus the extra en- 
tropy at constant time arising from the accreted 
matter; the last term denotes the deuterium lumi- 
nosity arising from both the protostar deuterium 
content and the one accreted with the infalling 
material. We adopt the same assumption as in 
Hartmann et al. (1997, see their Eqs(3)-(4)) and 
BCG09, namely that the accretion occurs over a 
fraction S of the stellar surface, so that 

L acc = 4nR 2 5F acc , (2) 

where F acc is the radiative flux averaged over the 
accreting region. Models of continuum emission 
from optical surveys suggest that S < 10% (Hart- 
mann et al. 1997), so we will assume S = in the 
rest of the paper, as in BCG09, for sake of sim- 
plicity. Further assuming, as in Hartmann et al. 
(1997), that there is a negligible energy transfer 
between the non accreting stellar surface and the 
region where accretion occurs, the stellar luminos- 
ity reads 

U = 4tt.R 2 (1 - <S)oT 4 ff w 4nR 2 aT^ (3) 

The impact of reducing the star radiating surface 
can be inferred from the study of Chabrier et al. 
(2007b), who explored the effect of magnetic spot 
coverage. These authors found that a surface spot 
coverage of 20% affects the radius evolution of 
young brown dwarfs by less than 5% (see Fig. 3 of 
Chabrier et al. 2007b). Consequently, for S < 10% 
or so, we do not expect a significant effect on the 
evolution/radius compared with the much larger 
impact due to accretion, as shown below. 

In Eq([T]) , the factor e is the amount of internal 
energy per unit mass brought up by the accret- 
ing material (e < 1 for spherical accretion while 
e < 1/2 if accretion occurs from a thin disk at the 
object's equator (Hartmann et al. 1997). As in 
BCG09, a value e = 1/2 will be assumed in the 
present calculations. Finally, the factor a < 1 is 
the fraction of accretion energy absorbed by the 
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protostar at the accretion shock, while (1 — a) is 
the fraction radiated away. The evolutionary mod- 
els are based on the same input physics and gen- 
eral framework as described in BCG09 and BC10, 
except for three main differences, as described in 
the next subsections. 

3.1. Accretion rates 

As mentioned above, the accretion rates are 
not arbitrary, but derived from the hydrodynamics 
simulations of prestellar core collapse mentioned 
in f|2j A supercritical shock at the protostar sur- 
face corresponds to a = 0, implying that in such a 
case all the accretion energy is radiated away and 
does not contribute to the protostar internal heat 
content. Such supercritical shock conditions have 
been found in ID complete RHD calculations of 
the first collapse (first Larson's core) of a 1 M Q 
and a 0.01 Mq prestellar core, both when using a 
grey approximation (Commergon et al. 2011) and 
a more accurate multigroup frequency description 
for the opacity (Vaytet et al. 2012). There is 
no garantee, however, that the same results still 
hold for the second collapse. Work to explore this 
issue is under progress. As shown in Hartmann 
et al. (1997) and BCG09, the value a crit ~ 0.1- 
0.2 corresponds to a critical value. While a value 
a < ctcrit leads to a contraction of the accreting 
protostar compared to the non-accreting counter- 
part at same mass and age, a value a > a cr u 
means that enough energy is absorbed by the ac- 
creting protostar to balance or even dominate its 
contraction, leading to a null effect or an expansion 
compared to the non-accreting counterpart. The 
first case refers to what is denominated " cold" ac- 
cretion, meaning that the entropy of the infalling 
material is significantly lower then the one of the 
protostar, while the second case corresponds to 
"hot" accretion, as briefly mentioned in SjTJ 

3.2. Initial protostar mass 

In BCG09, the protostar initial masses, that 
will be denoted M\, were chosen to range from 
0.001 to 0.1 Mq. This range was meant to bracket 
the typical value of the first Larson's core, as found 
in ID (Larson 1969, Masunaga et al. 1998) and 3D 
(Commergon et al. 2010, Tomida et al. 2010) nu- 
merical simulations of the first collapse. A more 
relevant value for the initial mass of the proto- 
star, however, is the one typical of the second adi- 



abatic Larson's core, which forms when the cen- 
tral density reaches values around ~ 10~ 2 g cm~ 3 , 
characterising the so-called stellar density in cloud 
collapse calculations (Larson 1969, Masunaga & 
Inutsuka 2000). Note, however, that this mass re- 
mains ill-determined, as no 3D simulations of the 
second collapse including the appropriate detailed 
physics (EOS, radiation hydrodynamics, magnetic 
field) have been carried out yet. We thus must rely 
on existing calculations or order of magnitude es- 
timates. 

A crude estimate for the second core's initial 
mass can be inferred as follows. Assuming a poly- 
tropic EOS, P = K n p 1+1 ' n , the structure of the 
second core can be estimated as a scaled down ver- 
sion of the first core (Masunaga et al. 1998), i.e. 
from a second-order Lane-Emden equation. This 
yields for the central density and pressure: 



A n p 



3A n M 



Pc = W r 



GM 2 



(4) 



AttR 3 ' R 4 
where M and R denote the mass and the radius 
of the second core, respectively, and A n and W n 
are the appropriate polytropic constants (Chan- 
drasekhar 1957). Assuming a perfect gas EOS 
yields for the central temperature 

T _ |im fl P c _ AnG^mH W n M 
Pc 3fc B A n R ' 

where p, denotes the mean molecular weight, run 
the atomic mass unit and fee the Boltzmann con- 
stant. The second hydrostatic core will form when 
all molecular hydrogen is dissociated, i.e. when 
the temperature reaches about 5000 K (Saumon 
et al. 1995). This yields for the second core mini- 
mum mass: 

3fc B A n 



Mi 



M 
V M G > 



M m i n — 



^TrGpniH W n 
T 



R 



4.0 x l0- 4 (-^)- 1 (^^)( 

y 1.0' v 5OOOK n i? 



.65 x 10 



— 2 ( A* \ — l ' 



R 



y 1.0' V 5000K A A.U. ' 

(6) 



where the values have been estimated for a fully 
ionized monoatomic polytropic gas n = 3/2. 
Taking R « IRq for the second core's radius, 
as inferred from the location of the first accre- 
tion shock in numerical simulations of the sec- 
ond collapse (Larson 1969, Masunaga & Inutsuka 
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2000), eqn.([6]) yields about half a Jupiter mass 
for Mi. The most recent numerical determination 
of the minimum mass for opacity-limited frag- 
mentation by Boyd & Whitworth (2005) yields 3 
M, Iup . Whitworth & Stamatellos (2006) have de- 
rived a minimum mass for star formation, denoted 
as minimum mass for primary fragmentation, in 
the range 1-4 Mj up . 3D hydrodynamics or MHD 
barotropic numerical simulations of the collapse of 
molecular clouds to stellar densities (Bate 1998, 
Intsuka et al. 2010, Machida et al. 2010, 2011) 
yield a second core mass Mi ~ 1 Mj up . 

It is unclear, however, to which extent varia- 
tions in the initial cloud prestellar core conditions 
will affect this value. The thermodynamics of the 
collapsing gas, for instance, may have some im- 
pact on the final (protostar) core mass. For a 
polytrope, the internal entropy is related to the 
central pressure and density as Si nt oc ln(-fW) oc 

Pc 

h^iJM 1 / 3 ) oc InM 4 / 3 , where we have made use 
of eqn.(3) (e.g. Stahler 1988). Assuming a per- 
fectly adiabatic collapse, this implies a depen- 
dence of the second core mass upon the entropy 
Mi oc exp{|(Sj n t — So)}, where Sq denotes the en- 
tropy at some initial stage before the onset of the 
second collapse. Because of the thermal properties 
of the gas, however, the collapse is not expected 
to be perfectly adiabatic, as part of the contrac- 
tion work (PdV) will be evacuated by radiative 
cooling, implying dQ ^ 0. The gas may also get 
heated by the accretion luminosity onto the cen- 
tral object. Indeed, variations of the entropy at 
the first core stage are found in RHD collapse sim- 
ulations (Commercon et al. 2011). These varia- 
tions, however, remain modest, < 5%. Although 
the impact at the second core stage remains to 
be explored, this suggests that variations of the 
second core mass due to the collapsing gas ther- 
mal properties are likely to remain at the few per- 
cents level, except in particular situations where 
the gas has been heated significantly due e.g. to 
strong radiative feedback from a nearby massive 
star. Similarly, variations in the initial rotation of 
the cloud prestellar core are also found to barely 
affect the initial mass of the protostar (Machida, 
private communication), suggesting a limited im- 
pact of rotation on the second core mass. Mag- 
netic braking, however, is found to increase the 
limit of fragmentation and thus the minimum core 
mass during the first collapse (Commercon et al., 



2010). Even though, again, no calculation exists 
for the second collapse yet, the impact of magnetic 
field at this stage is likely to remain weak, as most 
of the flux is expected to get lost by ambipolar 
diffusion during the first collapse. In summary, it 
seems difficult to significantly modify the initial 
mass of the protostar, Mi, even though this issue 
remains to be examined by dedicated calculations 
of the second collapse. 

Given these results, and awaiting for more ac- 
curate determinations of the second core mass 
from 3D full RMHD simulations of the second 
collapse, we have elected in the present calcula- 
tions to take for the value of the initial protostar 
mass, at the very beginning of the accretion pro- 
cess, a value Mi = lMj up , with possible varia- 
tions up to 5Mj U p. It is essential at this stage to 
stress that HOK11 consider only one value for this 
seed mas^J namely Mi — 10Mj up . As mentioned 
above, such a value is representative of the typical 
first Larson's core, but is likely to be quite a too 
extreme value to examine the early evolution of 
low-mass stars and brown dwarfs under standard 
cloud conditions. 

3.3. From cold to hot accretion 

The third difference w.r.t BCG09 is that only 
two values for a ^= were considered in these cal- 
culations, namely a= 0.2 and a= 1 (see Table 
1 of BCG09). The value a= 0.2 was found to 
be the one for which the protostar's accretion in- 
duced contraction was balanced nearly exactly by 
internal heating induced expansion, defining the 
limit between "cold" and "hot" accretion. In the 
present study, we explore more values for a, be- 
tween and ~ 0.2, in order to quantify the impact 
of even a modest fraction of absorbed accreting en- 
ergy upon the protostar evolution. Note that such 
finite values for a (a / 0) are considered only for 
sufficiently high accretion rates, namely rates ex- 
ceeding some critical value M CI > lO _5 M0yr _1 , 
as inferred from Fu Ori events (see §1). 



3 Note that the definition of the initial mass used in our cal- 
culations and in the calculations of HOK11 is the same, 
representing the initial mass for the stellar evolution calcu- 
lations. 
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4. The HRD luminosity spread 

A significant spread in the HRD for low-mass 
objects has been observed in many young, ^Myr 
old, clusters and star forming regions. The exact 
amount of the spread is difficult to assess, as it can 
vary between the regions and is affected by various 
sources of uncertainties. The typical L-T c g range, 
however, corresponds to the region bracketed by 
the 1 Myr and ~ 10 Myr Baraffe et al. (1998) 
isochrones for non-accreting objects in the LMS 
and BD domain, although some observations sug- 
gest an even larger spread (see e.g. Bayo et al. 
2011 for the Lambda Orionis star forming region). 
This typical observational spread is illustrated in 
Fig. 1 of BCG09. In the present section, we will 
show that episodic accretion can lead to such a 
typical spread after 1 Myr. 

4.1. Lower part of the HRD 

Using the conditions described in £j3j we first 
explore the impact of accretion for objects in the 
cool part of the HRD, in the region T g < 3000 
K. HOK11 claim that, in this domain, it is not 
possible to produce a luminosity spread, except 
if one adopts unrealistically small values for the 
initial radius (thus entropy) of the protostar (see 
Figs. 6 and 7 of HOK11). 

4- 1.1. Variations of the initial protostellar mass 

Figure [T] shows sequences of evolution obtained 
from the collapse of an initial prestellar core M c 
=0.085 M , with P = 2.3x 10~ 2 , for various values 
of M; between 1 and 5 Mj up , for a — 0. The inset 
in the figure illustrates the accretion history of the 
collapsing parent cloud core. Note that, for such 
a small initial core mass, the accretion rate barely 
exceeds lO _6 M0yr _1 , as mentioned in with 
only two weak accretion bursts which occur in the 
early phase due to a high (yet within the obser- 
vational upper limit) value of /3. We stress again 
that the accretion history (bursts, rates) during 
the evolution of the protostar is the one obtained 
from the core collapse calculations, insuring com- 
plete consistency of the calculations. All evolu- 
tionary sequences terminate at ~ 1-2 Myr and the 
positions at 1 Myr are indicated by the solid cir- 
cles. At this age, the final object produced along 
these sequences is a brown dwarf with a final mass 
Mi =O.O65-O.O7M , depending on Mi. As seen 



in the figure, a range of initial protostar masses 
between 1 and 5 Mj up does lead to the sought 
typical luminosity spread after 1 Myr, for a = 0. 
As discussed in <|3j such a variation among the 
protostar initial masses at the second core stage 
can not be excluded at the present time. As seen 
in the figure, protostar initial masses as small as 
Mj ~ 1 Afj U p are required in order to produce 
objects as faint as the ones corresponding to the 
10 Myr non-accreting isochrone^ Since, assuming 
uniform core-to-star mass conversion due to out- 
flows, as suggested by observations (e.g. Andre 
et al. 2010), small protostar masses can only 
be produced by small initial cloud prestellar core 
masses, this in turn implies the existence of small 
(< 0.1 Mq) prestellar core masses. Such cores are 
starting to be discovered with Herschel (Andre 
et al. 2010, Konyves et al. 2010) and are pre- 
dicted by analytical gravo-turbulent theories of 
cloud fragmentation (Padoan & Nordlund 2002, 
Hennebelle & Chabrier, 2008, 2009). As seen in 
the figure, the obvious reason why HOK11 could 
not succeed producing a spread in the HRD for 
objects with final temperature T e g < 3000 K, i.e. 
Mr <0.1M Q , is that their initial protostar mass 
is set up to be M; =10 Mj up , which indeed leads 
to a final object close to the 1 Myr non-accreting 
isochrone. In order to verify this statement, we 
have performed a calculation for Mj =10 Mj up 
(for the same prestellar core mass M c = 0.085M Q ): 
we indeed obtain an object with T c g ~ 3000 K 
and a position at 1 Myr very close to the 1 Myr 
non-accreting isochrone, similar to the results of 
HOK11 (see black filled square in Fig. [TJ. This 
confirms the fact that we can reproduce the calcu- 
lations of HOK11, if we use the same initial condi- 
tions, and that the very reason for their disagree- 
ment with the BCG09 results is the use of one sin- 
gle value, Mi =10 Mj up , for the initial protostar 
seed mass, a very likely too large value for the typ- 
ical second Larson's core mass, representative of 
this seed mass, as mentioned earlier. An other in- 
teresting effect illustrated in the figure is that, for 
the present prestellar core mass M c =0.085 Mq, 
the impact of early accretion is significant even af- 
ter 1 Myr for initial masses Mj < 5Mj up , while the 
effect vanishes for larger initial masses after about 
this time. A quantitative explanation for such a 



4 Such faint objects are indeed observed in few Myr old clus- 
ters (see BCG09 and references to observations therein). 
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behaviour, which involves deuterium burning, is 
given in Appendix A. A comparison between the 
final mass of accreting vs non-accreting objects for 
a given (L, T c g) location in the HRD will be given 
in ^ 

Figure [I] illustrates the effect of the variation 
of the initial protostar mass M; on the evolution- 
ary sequence, for a given prestellar core mass. As 
shown in BCG09 and as just mentioned, for a 
given accreted mass AM, the smaller M- u thus 
the smaller the initial binding energy, the larger 
the impact of accretion on its structure (see expla- 
nation in Appendix A). These calculations show 
that, considering different values of M; between 1 
and 5 Mj up naturally produces the spread in the 
cool part of the HRD. 

We stress again that, as stated in BCG09, the 
impact of accretion on the structure of an ob- 
ject depends essentially on the total amount of ac- 
creted mass, AM, since the initial protostar very 
formation, and depends only weakly on the de- 
tailed accretion history, i.e constant or monotonic 
vs episodic accretion (see §3.2 of BCG09). This 
has been confirmed by HOK11. Besides the rea- 
sons mentioned in §1, the obvious reason in fa- 
vor of episodic accretion is that it takes a signifi- 
cantly shorter time for a star to build up its mass 
from an initial protostar mass through a few in- 
tense accretion bursts than from a steady, or more 
likely rapidly decreasing (see the inset in the fig- 
ure) "standard" accretion rate. 

Fig. [I] also illustrates the impact of a finite 
value of a, i.e. of some fraction of accreting en- 
ergy absorbed by the protostar, on its evolution. 



As mentioned in S3.3 we define M cr as some crit- 
ical accretion rate above which a is switched from 
(for M < M cr ) to some finite value a max (for 
M > M cr ). Based on the Fu Ori event men- 
tioned in §1, we arbitrarily define this critical rate 
as M cr = lCP 5 M yr -1 . As seen in the figure, 
for a protostar initial mass Mi = lMj up , a max 
~ 2 %, i.e. a tiny fraction of accreting energy 
absorbed by the protostar, is sufficient to coun- 
terbalance the contracting effect of mass accre- 
tion, yielding eventually a position in the HRD 
for the final object after 1 Myr similar to that of 
its non-accreting counterpart at same mass and 
age. Interestingly enough, a significantly larger 
value, a max = 0.2, leads basically to the same re- 
sult for such small initial masses. As explained 



in §3.1 of BCG09, this stems from the fact that, 
for a ^ (even for values as small as 0.02), the 
extra energy absorbed during the few bursts with 
M > M cr yields an evolution which proceeds at 
higher luminosity and larger radius than the non- 
accreting counterpart, which implies shorter ther- 
mal timescales tkh = GM 2 /(RL). After the last 
intense burst (at log t = 4.5 in Fig. [TJ, for in- 
stance, the thermal timescale of the accreting ob- 
jects with a ^ is shorter than the accretion 
timescale, i.e. tkh < TV, = M/M. This is in stark 
contrast with the case a = 0, for which tkh > t~m 
for the same burst intensities. As explained in 
Chabrier et al. (2007a) and BCG09, in the first 
case, the object quickly contracts and converges 
towards a location in the HRD very close to that 
of the non accreting counterpart. On the opposite, 
in the second case the thermal timescale is much 
longer, of the order of a few Myr, and only at this 
stage will the initial conditions due to early ac- 
cretion be forgotten. Such high accretion episodes 
leading to a ^ 0, however, should be rare in small- 
mass collapsing M c cores, at least according to 
our choice of M cr . Indeed, M hardly exceeds the 
value 10 _5 M Q yr _1 for such cores, even for large 
values of (3, as seen in the inset of Fig. [I] (see 
also VB10). Accretion bursts in low-mass mod- 
els occur only for cores with rather high values of 
rotational support, (3 > 2 x 10~ 2 , because high 
disk-to-star mass ratios are needed to trigger a 
burst. Note, however, that even a small number 
of such events is sufficient to significantly alter the 
radius, thus the luminosity of the accreting object. 
Had we adopted a slightly larger value of M cr (say 
3 x 10~ 5 A/ Q yr _1 ), our criterion for a non-zero 
value of a would never have been fulfilled for the 
M c model shown in Fig. [T] In that case, only a 
variation among the initial protostar masses, Mj, 
could lead to the typical luminosity spread. Note 
that, according to the BCG09 and present calcu- 
lations, large values of a>0.2 would imply the 
existence of very bright very-low-mass Fu Ori like 
events, located above and on the rightmost side 
of the 1 Myr non-accreting isochrone. The lack 
of such observations, so far, in the very low-mass 
domain of the HRD thus suggests that either they 
last a very short time, or they are intrinsically 
very rare, or even absent, indeed suggesting very 
small values of a in this (low-mass) domain. This 
diagnostic, if confirmed in the future, would thus 



8 



confirm on observational grounds the results ob- 
tained in the present and VB10 collapsing core 
calculations, namely that (i) the burst frequencies 
and intensities directly correlate with the initial 
core mass, and (ii) that hot accretion (defined as 
a>0.2) is unlikely to occur for low-mass cores, 
due to the lack, or at least the rarity of intense 
accretion bursts. 

4-1-2. Variations of the initial prestellar core 
mass 

In this subsection, we explore a different strat- 
egy. Given that the second Larson's core mass is 
essentially fixed by the entropy of the gas at the 
end of H2 dissociation, i.e. by the H2 tempera- 
ture dissociation and the density at which the gas 
becomes opaque (the aforementioned typical stel- 
lar density), and that this core mass should not 
depend much on the initial cloud conditions, the 
protostar initial mass will always have basically 
the same value. Fig. [2] illustrates in that case 
the evolutionary tracks obtained for the same ini- 
tial protostar mass, M{=1 Mj up , for accretion his- 
tories obtained from four different initial prestel- 
lar core masses, namely M c — 0.061, 0.075, 0.085 
and 0.1 Mq, for a=0. The four sequences pro- 
duce objects with final mass 0.03, 0.04, 0.065 
and 0.09 M Q respectively. As seen in the figure, 
the smallest M c produce objects at 1 Myr with 
a position in the HRD intermediate between the 
1 and 10 Myr isochrones, while the two largest 
M c yield objects with luminosities corresponding 
to the 10 Myr non-accreting isochrone. As ex- 
plained in Appendix A, these two quantitatively 
and even qualitatively different behaviors reflect 
the balance between the protostar's binding en- 
ergy and deuterium burning nuclear energy, for 
objects massive enough to ignite D-fusion. These 
results show that, for a fixed protostar initial mass 
M[=l Mj up , and within the effective temperature 
range T e g w 2500—3000 K, a variation of the cloud 
initial prestellar core mass can produce a luminos- 
ity and T e g spread, populating again the region in 
the HRD between the 1 and 10 Myr non-accreting 
isocrones. Needless to say, a combination of these 
results with the ones obtained in the previous sub- 
section for a given prestellar core mass and a varia- 
tion of the values of a will produce a nearly homo- 
geneous distribution of young, ~ 1 Myr old stars 
and BDs over the entire region between these two 
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Fig. 1. — Evolution of accreting sequences for ini- 
tial seed masses M\ varying from 1 to 10 Mj up 
and accretion rates derived from the collapse of 
a prestellar core of mass M c — 0.085M Q and 
P = 2.3 x 10~ 2 . Dashed-dotted lines correspond 
to a=0. Filled circles on each coloured curves give 
the position at 1 Myr. All tracks terminate at a 
few Myr. The black filled square indicates the po- 
sition at 1 Myr of a sequence starting with Mj=10 
Mj up and a=0. The triangles indicate the position 
at 1 Myr of models calculated with M\=l Mj up , 
a m ax=0.02 (filled triangle) and a max =0.2 (open 
triangle) (evolutionary tracks are not shown for 
these sequences for sake of clarity). The two long 
dashed (black) curves are the 1 Myr and 10 Myr 
isochrones of Baraffe et al. (1998) for non accret- 
ing models. The inset shows the accretion rate (in 
M Q yr _1 ) as a function of time (in yr). 



non-accreting isochrones. The case of BD's below 
0.04 Mq will be examined in §7. 

To summarize the present section, we can say 
that a luminosity spread in the lower, coolest part 
of the HRD (T cff ^ 3500 K, M <, O.1M ) similar to 
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Fig. 2. — Evolution of accreting sequences for 
four different initial prestellar core masses, namely 
M c = 0.061 (solid, magenta), M c = 0.075 (dash, 
blue), 0.085 (long-dash, red), and 0.1 (dash-dot, 
cyan) M , starting from Mj=l Mj up , with consis- 
tent respective accretion histories, for a=0. Filled 
circles give the position at 1 Myr. The two long 
dashed (black) curves are the same as in Fig[l] 



the observed one can be produced 

(i) either by starting with a spread of initial pro- 
tostellar masses in the 1-5 Mj up range, for a—0, 

(ii) or by starting with an unique initial pro- 
tostellar mass, M ; ~ lM Jup , but considering col- 
lapsing low-mass cloud prestellar cores of different 
masses - as indeed expected from the core mass 
distribution resulting from the parent cloud col- 
lapse - with some variation of a between and a 
few per cents. 

Or, obviously, by a combination of both sce- 
narios. Note, however that, in both cases, it is 
necessary to start with protostar initial masses as 
small as ~ 1 Mj up to obtain LMS or BDs as faint 



at 1 Myr as non-accreting 10 Myr old objects of 
same mass. The unambiguous observation of such 
faint low-mass objects in young (~ 1 Myr) clus- 
ters would thus indicate that the minimum mass 
for star formation indeed extends down to about 
this limit. In contrast, the demonstration that 
star formation can never produce such small initial 
protostar masses would suggest that either some 
physics is missing in the present calculations, or 
indeed the luminosity spread in the low mass part 
of the HRD is not, or at least not entirely due to 
early accretion in the protostar phase, or that the 
observations are questionable. 

4.2. Upper part of the HRD 

We now focus on the upper, hotter part of 
the HRD, namely final objects with T e g > 3000 
K, up to ~ 1 M Q . Always assuming a roughly 
mass-independent fraction of core accreting en- 
velope removed by outflows, such large masses 
must be issued from larger cloud initial prestellar 
core masses. Figure [3] displays sequences of evo- 
lution obtained for an initial prestellar core M c 
=1.53 M & , with /3 = 1.3 x 10~ 2 , starting from 
a protostar initial mass M\ = lMj up , for differ- 
ent values of a max (we recall that a is switched 
from to a max only during burst events exceeding 
M > M a = lO _5 M yr- 1 ). The prestellar core 
collapse calculations were performed up to 4=1.5 
Myr. At this age, the star has a mass M=0.75 M 
and still accretes at a rate M~8x lO _8 M0yr _1 . 
At 1 Myr, the central object has a mass M=0.71 
Mq. The results show that variations of a max in 
the range 0-0.2 produce the typical observed lumi- 
nosity spread in the HRD. Note in passing the in- 
credibly erratic evolutionary path for a max = 0.2. 
These tracks, as well as the ones illustrated in the 
previous figures, highlight the complex early evo- 
lution of accreting protostar until they reach their 
final mass, quite different from the usual const ant- 
T c ff Hayashi track typical of low-mass non accret- 
ing objects. For illustration, we also display the 
position at 1 Myr of a sequence starting with an 
initial mass M\ = 5 Mj up and a max =0 (empty cir- 
cle) . As seen, starting with M\ = 1 or 5 Mj up 
leads basically to the same results, a consequence 
of the large amount of total mass accreted and the 
resulting large binding energy of the protostar. In- 
deed, for such sequences, the energy release from 
deuterium burning can never overcome the gravi- 
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tational energy increase (as explained in Appendix 
A), and the effect of accretion is the same, inde- 
pendently of the initial mass. A similar spread 
was obtained, with cold accretion, in BCG09, for 
larger values of Mj (see discussion in S 4.1.1 ). As il- 



lustrated in Figure [3j adopting values for Mi more 
consistent with the ones expected for the second 
Larson's core, the only possibility to counterbal- 
ance the contracting effect of mass accretion on 
the protostar and to produce a luminosity spread 
is to assume a spread in the values of a between 
and ~ 0.2, i.e. some variations in the shock con- 
ditions on the protostar. This is consistent with 
the accretion histories obtained in the present and 
VB10 calculations, as well as in 3D simulations 
(Machida et al. 2010). Indeed, in contrast to the 
case with small prestellar cores, the accretion rate 
for the larger cores (even with moderate values 
of j3 > 5 x 10~ 3 ) often exceeds our critical value 
M cr = 10~ 5 A/ Q yr _1 and even 10" 4 Meyr^ 1 dur- 
ing the most intense bursts. It is intuitively ex- 
pected that larger accretion rates, thus larger ram 
pressures and hotter accreting material, are more 
likely to affect the internal heat content of the ini- 
tial protostar, fixed by the second Larson's core 
adiabat. In Appendix B, we provide some ana- 
lytical estimates of the properties of the accretion 
shock as a function of the accretion rate which 
provide some physical foundation for such larger 
values of a for higher mass objects. 

5. Fu Ori bursts 

As displayed in the inset of Fig. |3j the calcu- 
lated accretion rates corresponding to the collapse 
of a prestellar core of mass M c =1.53 M exhibit 
strong variations with time, with violent bursts 
reaching values characteristic of Fu Ori-like bursts. 
As mentioned above, such high rates are never 
reached during the collapse of the smaller cores 
investigated in |4.1| In case of strong accretion, 
approaching or exceeding the hot accretion limit, 
the most intense bursts strongly and durably affect 
the structure of the protostar, provoking its rapid 
expansion followed by a rapid contraction as soon 
as the burst stops. This is illustrated in Fig. |1J 
which portrays the evolution of the radius for pro- 
tostars of same initial mass M\ — lMj up and final 
mass M{ = 0.75M Q , accreting the same amount of 
mass AM, but for different values of a max , namely 
0, 0.04 and 0.2. The most intense accretion burst 
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Fig. 3. — Evolution of accreting sequences with Mi 
=1 Mj up and different values of a, with accretion 
rates (plotted in the inset in M Q yr _1 versus time 
in yr) derived from the collapse of a prestellar core 
of mass M c = 1.53M Q and /? = 1.3 x 10~ 2 . Filled 
circles on each coloured curves give the position 
at 1 Myr. The (magenta) open circle close to the 
filled circle corresponds to the position at 1 Myr 
of a sequence starting with M; = 5 Mj up and with 
amax^O. The two long dashed (black) curves are 
the same as in Fig. [T] 



for this sequence reaches 7.5 x 10 _4 Af Q yr _1 at an 
age of 0.12 Myr, when the protostar has built up 
a mass of ~ 0.35 M©. Despite the extremely short 
duration of the burst (~ 20 yr) , typical of Fu Ori 
bursts, absorption of 20% of the accreting energy 
by the protostar yields an expansion of up to a 
factor 2 in radius (blue dotted curve in Fig. [4j, 
reaching ~ 4 R Q . Most interestingly, this sequence 
matches remarkably well the properties of Fu Ori 
derived by HZC11, namely a central star mass ~ 
0.3 M Q , an inner disk radius ~ 5 R Q and an accre- 



tion rate M - 2 x 10 _4 M 



©yr 



Our calculation 
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thus explains the factor ~ 2 increase in the stellar 
radius inferred from the spectral energy distribu- 
tion of Fu Ori objects. As mentioned in Sj2] such 
high burst intensities require massive disks, very 
prone to strong and frequent episodes of gravi- 
tational instability and fragmentation, and thus 
massive (large M c ) prestellar cloud cores. Such 
cores do not need to rotate extremely fast and 
moderate values of j3 > 5 x 10~ 3 are sufficient 
to produce multiple bursts. However, further de- 
creasing /3 would result in disks of too small size 
(due to the corresponding decrease in the centrifu- 
gal radius) that are unable to fragment and such 
models would stop producing bursts. This is in- 
deed what is found in the hydrodynamics collapse 
calculations (sec e.g. VB10 and V10). 
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6. The global scenario 

The results presented in the previous sections 
suggest the following global scenario for the early 
evolution of accreting protostars and proto brown 
dwarfs. Objects in the LMS/BD mass regime 
form from the collapse of prestellar cores of various 
masses M c , as produced by the so-called core mass 
function (CMF), and various initial rotational to 
gravitational energy ratios /3 (and magnetic flux 
intensity as well). There is indeed ample evidence, 
in particular with the recent results of the Herschel 
mission, that the bulk of the stellar initial mass 
function (IMF) resembles the CMF and thus that 
most stars/BDs should form from the collapse of 
their parent core (Andre et al. 2010, Konyves et al. 
2010), in agreement with predictions of analytical 
theories of gravoturbulent fragmentation (Padoan 
& Nordlund 2002, Hennebelle & Chabrier 2008, 
2009, Chabrier & Hennebelle 2010, 2011). Initial 
protostar seed masses, resulting from the second 
collapse, are assumed to be of the order of M\ ~ 1 
M]up, as expected from analytical estimates and 
as found in existing numerical simulations. Some 
variations in Mi, however, depending on the ini- 
tial cloud conditions (e.g. magnetic field intensity 
or orientation) are presently not excluded. Fur- 
ther studies are definitely needed to explore this 
issue. The larger and the more rapidly rotating 
the initial cloud core, the more massive the cen- 
trifugal disk mass around the central protostar, 
and the more intense and frequent the accretion 
burst episodes. As shown above, the combination 
of high M c and/or large j3 can produce Fu Ori 
type bursts (M > 10~ 4 M Q yr -1 ) and explain Fu- 
Ori properties (accretion rate, mass, radius). A 
prediction of such a scenario is that there should 
be a dearth of Fu-Ori like systems for very low- 
mass objects (i.e mostly for masses in the BD 
regime). Assuming, as intuitively expected and as 
supported by the analytical arguments presented 
in Appendix B, that there is a correlation between 
a, the fraction of accreting energy absorbed by 
the protostar, and the burst intensities, we predict 
the accretion process onto the central object to be 
"cold-like" (a=0 or at most a few percents) for 
low accretion rates {M < lO _5 Af yr _1 or so), i.e 
for prestellar cores with small M c or small /3, and 
to become "milder" (0 < a < 0.2) with increas- 
ing M c and /?, possibly experiencing eventually 
some episodes of intense "hot" accretion bursts 
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(a > 0.2), characteristic of Fu-Ori like events. 
Such a link between a and accretion rates remains 
to be demonstrated on robust grounds, and will be 
briefly addressed in S|8] 

As mentioned above, for low mass prestellar 
cores, for which burst events do not occur and 
the accretion rate remains weak (<C lO -5 M0yr -1 ) 
most of the time, the accretion process is predicted 
to always be essentially " cold-like" . The spread of 
luminosities after ~ 1 Myr in the cool part of the 
HRD, resulting in an apparent age spread, can be 
explained either by variations of the protostar ini- 
tial mass within the range 1 Mj up < M\ < 5 Mj up , 



or, for a fixed initial mass M, 



1M 



Jupi 



by a 



distribution of initial prestellar core masses and 
some modest variations of a (0 < a < 2%). For 
more massive prestellar cores, the sought luminos- 
ity spread is best explained, for fixed protostar 
initial mass, by a larger spread in a values, up 
to 20% or so, as mentioned above. It is worth 
mentioning, at this stage, that there is not nec- 
essarily a direct correlation between the intensity 
of the bursts and the degree of luminosity spread. 
As explained in Sect. 4. 1.1 and appendix A (see 
also BCG09), the impact of accretion depends es- 
sentially on the fact that the accretion timescale, 
Tjv = M/M , is shorter or not than the thermal 
timescale, txH- As shown in the previous section, 
very intense bursts may lead to situations such 
that txH < t m with large values of a (not men- 
tioning the role of D-burning) , yielding very small 
or no luminosity spread, whereas more moderate 
accretion rates could lead to the opposite situa- 
tion. It also depends on the initial core mass, M c 
(see appendix A). 

This global scenario leads to a significant lu- 
minosity spread over the entire low mass part 
(< 1 ill©) of the HRD at 1 Myr, as illustrated for a 
few cases in Fig. [5]for various initial prestellar core 
conditions (M c , f3), protostar seed masses (Mj) 
and absorbed fraction of accreting energy (a). It 
is obvious that variations among these parameters 
will populate the entire HRD regions between the 
~1 and ~10 Myr non-accreting isochrones in the 
low mass star/BD regime. Work is in progress to 
perform a more systematic study, assuming a dis- 
tribution of M c , f3 and M\. As mentioned above, 
the collapse of cloud cores with large M c and /3 
produce the most intense accretion bursts, which 
can explain the observed properties of Fu Ori like 



objects, in particular their SED and the factor 
2 increase of the stellar radius. 



-2 



^ \ □ 

0.70 Mi* ^\ 

% X 

0.50 M N 



\ 



\ 



\ \ 
0.26 U m # § \ 

\ ^ 
0.09 M \^ 

0.065 



\ 



0.04 M. \ »\ 

\ \ 
0.03 M ^ x v 

X 





4500 



4000 



3500 3000 



2500 



Fig. 5. — Illustration of the luminosity spread in 
the HRD obtained for various sequences of evo- 
lution calculated for prestellar core masses M c = 
0.061, 0.075, 0.085, 0.1, 0.38, 0.92, 1.53 and 1.85 
Mq, leading to objects in the mass range 0.03 - 
0.9 Mq at 1 Myr, as indicated in blue on the left 
hand side of the symbols. The symbols give the 
position at 1 Myr. Solid (blue) symbols corre- 
spond to Mi=l Mj up ; open (red) symbols corre- 
spond to M[=5 Mj up . Circles: a max =0; squares: 
amax=0-02; triangles: a max =0.2. Black dashed 
curves are the same as in Fig. [T] 



7. Impact on the final mass and IMF de- 
termination 

A natural question arising from the present 
global scenario is : what is the impact of early 
accretion upon the mass determination at young 
ages and thus what is the impact on the IMF 
determination for young clusters, derived from 
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mass-luminosity relations based on non-accreting 
isochrones ? 

A first interesting feature concerning the im- 
pact of accretion upon the final object is that 
this impact is relatively small in the low-mass 
(<0.04Af Q ) BD domain. First of all, as seen in 
all figures, the dispersion of luminosity in this do- 
main reduces significantly, as the 1 and 10 Myr 
non-accreting isochrones tend to merge towards a 
single line. This stems from the fact that the lumi- 
nosity barely evolves during this timescale for low- 
mass BDs with masses <0.04M Q , due to ongoing 
deuterium burning] A second (related) feature, 
as seen in Fig. [5j is that accretion has a lesser 
impact on the evolution of BDs of final masses < 
0.04 Mq. As explained in Appendix A, the en- 
ergy released from (initially present and accreted) 
deuterium fusion can partly overcome the gravi- 
tational energy increase due to mass accretion for 
these very low masses, limiting the total contrac- 
tion of the object compared to its non-accreting 
counterpart of same mass and age. We expect this 
to be true as long as deuterium fusion remains ef- 
ficient, i.e down to ~ 0.01 Mq. Therefore, brown 
dwarfs in the mass range ~ 0.01 — 0.04Af Q should 
be only moderately affected by their accretion his- 
tories. Work is under progress to explore the early 
evolution of accreting BD's below the deuterium 
burning minimum mass. 

This is different for higher masses. Table 2 of 
BCG09 already provides some information about 
the difference between observed properties, L and 
T e g, for accreting and non-accreting objects of 
given mass in the LMS domain. In order to il- 
lustrate the impact of accretion history on the 
IMF determination, we provide some comparisons 
in Table 1 and Table 2 between the masses in- 
ferred from accreting and non-accreting (Baraffe 
et al. 1998) models, respectively, for a few typical 
examples. Table 1 compares the inferred masses 
for given luminosity L and effective temperature 
T e s while Table 2 compares the masses for given 
age and effective temperature T c g. As seen in Ta- 
ble 1, inferring the mass for a given (L, T e g) from 
non-accreting models can yield an overestimation 
of the mass by ~ 20% to ~ 40% for the largest 



5 Dcutorium burning lasts ~ 9 Myr for a 0.03 Mq BD and 
~ 6 Myr for a 0.04 Mq BD according to the Baraffe et al. 
(1998) models. 



masses. The overestimation can be even larger 
when inferring the mass from given age and T e g, as 
illustrated in Table 2. This is consequential when 
inferring the IMF for young clusters from the ob- 
served signatures in the HRD. Similarly, as seen 
in the same Table, inferring an age for LMS from 
their observational signatures from non- accreting 
isochrones can lead to ages as large as 20 Myr, 
whereas the objects are 1 Myr old in reality. This 
highlights the crucial impact of early accretion his- 
tory on the evolutionary properties of young low- 
mass objects. As mentioned above, for masses 
M^O.O4M , the impact is much less severe. A 
more systematic study exploring a wider range of 
parameters (M c , /3, M\, a) is under progress to 
further explore this impact on the IMF. 

8. Conclusion 

In this paper, we have presented a global sce- 
nario, based on accretion evolutionary sequences 
for low-mass stars and brown dwarfs characterized 
by accretion histories derived from numerical hy- 
drodynamics simulations of their parent collapsing 
cloud prestellar core. This scenario suggests the 
following paradigm: 

• (i) cloud fragmentation and collapse lead to 
a distribution of collapsing prestellar cores 
with various masses M c , various ratios of ro- 
tational to gravitational energy /3, and vari- 
ous magnetic field properties. Each of these 
collapsing cores leads to different episodic 
accretion histories, characterized by differ- 
ent burst episodes/intensities. As a general 
feature, the larger M c and /?, the more prone 
the disk to fragmentation and the more in- 
tense the accretion burst episodes. 

• (ii) initial protostars (second Larson cores) 
form at the center of these collapsing prestel- 
lar cores when reaching the stellar density 
(~ 10~ 2 g cm~ 3 ) and thus the adiabat de- 
termined by H2 dissociation and the opacity 
limit for fragmentation. This limit remains 
presently uncertain but should be of the or- 
der of 1 Mj U p and should vary only mod- 
estly with the cloud conditions, even though 
the exact amount of variation remains to be 
determined from appropriate RMHD simu- 
lations of the second collapse. 
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Table 1: Comparison between masses (in Mq) inferred from accreting and non accreting models, respectively, 
for some given (L, T c g). The table provides the mass MiMyr resulting from accreting sequences which 
correspond to these (L, T e ff) at 1 Myr, the mass Mbcah98 which corresponds to the same (L, T e g) for non 
accreting models (Baraffe et al. 1998), with an age iBCAH98 ( m Myr). The last three columns provide the 
various initial conditions, M c , Mi and a, which yield these masses for the accreting sequences. 



log L/L Q 


T c tf 


Ml Myr 


MbCAH98 


£BCAH98 


M c (Mq) 


M (M Jup ) 


a 


-2.12 


2783 


0.04 


0.04 


5 


0.075 


1 





-1.919 


3026 


0.065 


0.08 


10 


0.085 


1 





-1.72 


3115 


0.09 


0.11 


10 


0.1 


1 





-0.83 


3750 


0.5 


0.7 


20 


0.92 


5 





-0.74 


3763 


0.5 


0.7 


15 


0.92 


1 





-0.32 


4182 


0.9 


1.1 


10 


1.85 


1 





-0.37 


4221 


0.9 


1.05 


15 


1.85 


5 






Table 2: Same as Table [T] for a given (age, T c g). M acc is the mass derived from an accreting sequence with 
M c , Mj and a specified in the table. 



age (Myr) 


T e ff 


Macc 


MbcAH98 


M c (M ) 


M (M Jup ) 


a 


1 


2720 


0.03 


003 


0.061 


1 





1 


3026 


0.065 


0.11 


0.085 


1 





1 


3750 


0.5 


0.9 


0.92 


5 





1 


4182 


0.9 


1.4 


1.85 


1 






• (iii) a fairly small or negligible (at most a few 
%) amount of accretion energy is expected 
to be absorbed by the central object in the 
very-low mass (mostly BD) regime, as ex- 
pected from the moderate accretion rates 
produced by small-mass collapsing cores, 
which yield very scarce of even no burst 
episodes. In contrast, a larger amount of 
absorbed accretion energy (up to ~ 20%) 
is expected for larger initial core masses 
(thus larger stellar masses) and large rota- 
tion rates. This amount increases with the 
number of accretion bursts with rates ex- 
ceeding some critical value ~ 10~ 5 Mq yr _1 . 
Substantial variations of a in this regime are 
expected, with a increasing with the core 
mass and initial rotation rate. 

• (iv) in some cases, high-mass and rapidly 
rotating collapsing cores may lead to such 
intense bursts that one enters the domain 
of "hot accretion", with a > 0.2, yield- 
ing a strong increase of the protostar radius 
and luminosity, typical of Fu Ori like events. 



Whether the evolution of a from "cold" to 
"hot" is determined by a change in shock 
conditions or by a change of topology in the 
accretion process remains to be determined 
and should motivate dedicated studies. 

We have shown that this global scenario can ex- 
plain both (1) the observed luminosity spread in 
the HRD of young clusters over the entire low- 
mass (< 1M ) range, with no significant age 
spread, and (2) the properties (mass, radius) in- 
ferred from the spectral energy distribution of 
Fu Ori objects, providing an explanation to the 
questions raised by HZC11. If correct, this uni- 
fied picture should close the controversy raised 
by HOK11, except if it can be shown unambigu- 
ously that the initial protostar/BD mass can not 
be smaller than about 10 Mj up , i.e. that the first 
and second Larson's core are about the same mass, 
a rather unlikely possibility given the enormous 
differences in temperatures and densities, thus en- 
tropies, during the first and second collapse. 

More work is definitely required to test the 
above-mentioned paradigm. Regarding assump- 
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tion (i), the link between the accretion burst in- 
tensity and both M c and j3 requires further stud- 
ies, ideally with 3D MHD simulations including 
the effect of magneto-rotational instability in the 
innermost disk, but as a first step by including 
self-consistent radiative feedback from the proto- 
star as discussed in section [2] Concerning point 
(ii) , 3D simulations of the second collapse includ- 
ing radiation hydrodynamics and non-ideal MHD 
are now underway (Machida et al. 2010, Masson 
et al. 2011). They should provide more insight 
on the second Larson's core mass determination 
and characteristic mechanical and thermal prop- 
erties (e.g initial entropy content) and on the de- 
pendence of these properties upon the cloud initial 
conditions (e. g rotation, magnetic field). 

Assumptions (iii) and (iv) suggest the existence 
of a threshold in accretion rates above which a 
transition from cold to hot accretion should occur. 
Many possibilities can be speculated to suggest an 
evolution of the accretion pattern/properties onto 
the protostar. The simplest explanation, as out- 
lined in Appendix B, is that large accretion rates, 
orders of magnitude larger than the "standard" 
c|/C? ~ 10~ 6 M©yr _1 value, are energetic enough 
(1) to penetrate substantially below the surface 
layers, within the protostar deeper (convective) 
layers, (2) to dissipate enough kinetic energy to 
significantly affect the protostar's heat budget, 
halting its contraction. Such a transition may 
also stem from a change in the accretion geometry, 
evolving for instance from magnetospheric accre- 
tion for weak accretion rates, to thick disk accre- 
tion for intense bursts, as suggested by the absence 
of magnetospheric accretion in Fu Ori (HZC11). It 
might also stem from a change in the disk proper- 
ties, evolving from thin disk to thick disk accretion 
with increasing accretion rates. It could also cor- 
respond to a transition from a supercritical shock 
to a (optically thin) subcritical shock at the sec- 
ond core stage (Commeroon et al. 2011). The 
topology (and the value) of the protostar's mag- 
netic field may also play a major role in the accre- 
tion process; 3D MHD simulations by Long et al. 
(2008) indeed show that the mass accretion rate 
from the disk onto the star, the area covered by 
hot spots and the angular momentum flux from ac- 
creting material between the disk and the star are 
several times smaller in the case of a quadrupole 
dominated configuration than in the dipole case. 



The dependence of the protostar magnetic field 
topology upon its mass could thus have an impact 
on the value of a. Finally, increasing the accre- 
tion rate yields a smaller magnetospheric radius, 
at which the disk is truncated by the stars's mag- 
netic field, evolving from accretion through funnels 
to direct accretion onto the star through a bound- 
ary layer (Romanova et al. 2008a). It is indeed 
known that at high accretion rates characteristic 
of Fu Ori outburst phases, the magnetospheric ra- 
dius moves much closer to the stellar surface. High 
accretion rates are also predicted to yield unsta- 
ble situations with accretion evolving from funnel 
streams to a situation with multiple hot spots on 
the stars's surface (Romanova et al. 2008b, Kulka- 
rni & Romanova 2008). Exploring these issues re- 
quires dedicated observational and theoretical ef- 
forts in order to better understand the star-disk 
interaction and the accretion processes onto pro- 
tostars. 

The consequences of the present study have 
also far reaching consequences in the domain of 
IMF and star formation. First of all, it shows 
that the concept of a birthline does not apply to 
low-mass (< 1M Q ) objects. Each proto-star/BD 
experiences its own (unpredictable) accretion his- 
tory and ends up randomly in the HRD at the 
end of the accretion process. Second of all, infer- 
ring masses from observed L-T e g or T c g-age prop- 
erties in the HRD from isochrones which do not 
take into account the complete accretion history 
can yield severely incorrect determinations, possi- 
bly overestimating the mass by as much as 40% 
or more. The possibility of such errors should be 
borne in mind when trying to determine the IMF 
of young (< 10 Myr) clusters. In the same vein, 
age determinations of faint objects in clusters or 
star forming regions from "standard" isochrones 
can overestimate the age of these objects by sev- 
eral tens of Myr ! As we have shown in BC10 
that accretion also strongly increases lithium de- 
pletion, this definitely invalidates the arguments 
in favor of "slow" star formation and should close 
this debate. 

To conclude, we are conscious that the scenario 
presented in this paper (and in BCG09 and BC10) 
has weaknesses, as it relies on some assumptions 
which, so far, lack robust physical fondations. On 
the other hand, the scenario relics on the treat- 
ment of the impact of episodic accretion, inferred 
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from the collapse of a cloud prestellar core, onto 
the structure and evolution of the accreting pro- 
tostar/BD. Within this framework, it provides 
a consistent explanation for several, apparently 
uncorrelated, observational features, namely the 
HRD luminosity and radius (see Jeffries 2007) 
spread, the properties of Fu Ori, the abnormal 
Li depletion of some young objects (see BC10) 
and the luminosity distribution of embedded pro- 
tostars (Evans et al. 2009, Dunham & Vorobyov 
2011). This should motivate further observational 
and theoretical studies to test the reliability of this 
global picture. 
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A. Mass and deuterium accretion 

In this appendix, we explain the results displayed in Fig. [T] namely the increasing effect of accretion on the 
central object's structure for decreasing Mi at fixed M c , for M c = O.O85M . These results can be understood 
from arguments suggested by Lee Hartmann, by comparing the binding energy of the central object with the 
amount of energy available from deuterium fusion (Hartmann, priv. com.). Assuming n = 3/2 polytrope, 
the total binding energy of a protostar is: 

E ^ = -7^r- (A1) 

The fusion reaction of deuterium p + d provides an energy Q—5A9A Mev per reaction (Caughlan & Fowler 
1988). Assuming a mass fraction [D] = 2 x 10~ 5 , the energy per gram of material available from D burning 
is thus edeut ~ 10 14 erg/g. Assuming complete fusion of D (initially present and accreted), the total energy 
released from D fusion when the object has reached a mass M, is: 

-Edcut = M x e dout . (A2) 

If |i?tot| > -Edeut, the release of nuclear energy is insufficient to balance the increase of gravitational energy 
of the growing protostar, yielding a more compact object than its non-accreting counterpart of same mass 



M and age (Stahler 1988; Hartmann et al. 1997). Combining Eqs (Al) and (A2), the above-mentioned 
inequality yields the following condition on the radius of the protostar: 

M 

R < 8.2 — Rq = R clit . (A3) 

Mq 

If R < i?crit, accretion produces an object more compact than the non-accreting counterpart, provided 
t m < Tkh- The latter condition is always fulfilled during burst events (with a=0). In contrast, if R > -R C rit, 
the energy released by D burning will limit the contraction of the accreting object. For the models displayed 
in Fig. [lj comparison of the models starting with 1 Mj up (hereafter model A) and 5 Mj up (model B), 
respectively, shows that before the onset of D fusion, model A contracts faster than model B as a result of 
accretion, for the same accretion rate, on an object with an initially lower binding energjj^] and no nuclear 
energy generation. At the onset of deuterium burning (T c ~ 6 x 10 5 K, t = 10 4 yr), model A has a larger 
binding energy than model B, because of its significantly smaller radius (see Table[3]). Deuterium burning and 
ongoing accretion thus proceed with R < R CI it for model A and with R > i?, cr it for model B. Consequently, 
according to the afore mentioned argument, accretion has a much larger impact on model A than on model 
B, yielding a significantly more compact structure than the non-accreting counterpart (see Table [3]). 



The condition ( A3 ) thus provides a good criterion to infer the effect of accretion during and at the end 
of the accretion process, for objects massive enough to ignite D-fusion, i.e. M^0.013M Q (Chabrier et al. 
2000). Note that for model A, at 1 Myr the object still has a radius R < R CI n, hence looking older than the 
non-accreting counterpart, whereas for model B the object has a radius R > R C nt at this epoch, closer to 
that of the non-accreting counterpart (see Table |3j. 



Criterion (A3) explains as well the results displayed in Fig. [2] (see [4.1.21, namely the larger the prestellar 



core mass M c the more compact the structure of the newly formed object compared to its non-accreting 
counterpart at 1 Myr (see Table Eh. Indeed, for a given (same) initial protostar mass M;, larger prestellar 
core masses M c yield objects with a larger final mass and thus a larger binding energy, so that criterion 
(A3) is more easily fulfilled (see Table EJ). As mentioned above, these arguments hold as long as the forming 



object is massive enough to fuse deuterium. 



6 Both sequences start with an initial radius R ~ IRq. Variations of the initial radius for values > 0.7 Rq yield the same 
conclusions. 
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Table 3: Properties of accreting models displayed in Fig. [I] (see (4.1.1) and starting respectively with 1 Mj up 



(model A) and 5 Afj up (model B), for a—0. Models with M = refer to the non-accreting counterparts 



Energies are in units of 10 45 erg. The Kelvin-Hclmholtz timcscale tkh and accretion timescale tv, are also 



provided (for the latter, only for models A and B at 10 yr when accretion is relevant), in units of Myr. 



model 


t (yr) 


M (M ) 


R(R G ) 


-Rcrit (-R©) 


Etot 


^tdout 


tkh (Myr) 


(Myr) 


A 


10 4 


0.036 


0.3 


0.3 


7.2 


7.2 


25 


0.07 


B 


10 4 


0.04 


0.52 


0.33 


5.1 


8.1 


6 


0.08 


A 


10 6 


0.065 


0.4 


0.54 


17.5 


13.1 


27 




B 


10 6 


0.07 


0.65 


0.57 


12.2 


13.9 


8.5 




M = 


10 6 


0.065 


0.70 








6 




M = 


10 6 


0.07 


0.75 








6 





core masses M c and same Mj 


= lM Jup 


, and of non-accreting counterparts. 




model 


M (M ) 


R (R Q ) 


-Rcrit {Rq) E to t 


(10 45 erg) £ tdcut 


(10 45 erg) 


M c = O.O61M 


0.03 


0.35 


0.26 


4.84 


6.46 


M c = 0.075M© 


0.04 


0.38 


0.33 


7.1 


8.1 


M c = 0.085M© 


0.065 


0.4 


0.54 


17.5 


13.1 


M c = O.1M 


0.09 


0.48 


0.74 


28 


18 


M = 


0.03 


0.49 








M = 


0.04 


0.5 








M = 


0.065 


0.7 








M = 


0.09 


0.91 









B. Properties at the shock front 

In this appendix, we derive analytical estimates for the various properties of the flow at the shock front. 
All quantities are expressed in cgs units. 

• Pressure of the accretion flow 

The ram pressure of the accretion flow, falling down on the proto-star/BD at the free-fall/Keplerian 
velocity u acc = (2GM/i?) 1 / 2 , can be estimated as : 

_ M v &cc ^ 3 x 10 4 M M , R /2 

acc ~ 8nR?8 6 l lO- 6 Af yr-i n Af J l i? j 1 ' 

Assuming a stellar surface fraction S w 0.1 covered by accretion, this yields P acc ~ 10 4 ( 10 -e M 0yr -i ) dyne 
cm~ 2 for a 1 to 10 Afj up initial protostar. 

• Properties of the photosphere 

For an atmosphere in hydrostatic equilibrium, VP = — pg, where g — 2.7 x 10 4 (]^)(j^)~ 2 is the surface 
gravity, the optical depth at depth z, t(z) — z/l U) where l v = l/(pR) is the photon mean free path and R 
the Rosseland mean opacity, obeys the condition 

t(z)= / -dP. (B2) 
Jz 9 

The pressure and density at the photosphere, r ~ 1, are thus given by 
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where p denotes the mean molecular weight and 1Z the perfect gas constant. At ~ 10~ 6 Af yr _1 , the lMj up 
proto-star/BD core quickly reaches a few Mj up , which yields T c r > 2000 K and, for an initial radius ~ IP©, 
logg ~ 2.0-2.5 from ^a few to about 10 Mj up . In this regime, for R — p/T§ ~ 10-100, where Tq is the 
temperature in units of 10 6 K, the Rosseland mean opacity is dominated by molecular lines and is R sa 10~ 2 
cm 2 g" x (Fergusson et al. 2005). This yields P p h ~ 10 4 dyne cm -2 for the aforementioned mass range, in 
good agreement with the values obtained with the complete calculations (see Fig. 4 of Baraffe et al. 2002). 
• Low accreting rates 

Let first consider small collapsing prestcllar cores (M c < 0.1 M©). Assuming a ~30%-70% core-star mass 
conversion efficiency, these cores will produce final objects in the BD domain. As mentioned in the text, these 
collapsing cores do not produce high accretion bursts and the accretion rate remains essentially of the order 
of 10 _6 M©yr _1 (see e.g. inset in Fig. 1). According to the above estimates, the pressure at the accretion 
shock thus remains of the order of or at most about an order of magnitude larger than the photosphcric 
pressure, as the object builds up its mass, from ~ a few Mj up to, eventually, a few tens of Afj up . As will be 



shown below (see eqn.(B5l), for such conditions, the shock is likely to remain nearly superficial and thus to 
entail a negligible amount of energy compared with the proto-star/BD binding energy. This leads support 
to a "cold accretion" condition, with a value a <C 1 for the fraction of absorbed accreting energy in the very 
low mass domain. 

• High accreting bursts 

As the mass of the collapsing cores, and thus of the final objects, increases, the simulations predict an 
increasing number of bursts, many of them reaching the 10 _5 -10 -4 M©yr _1 level, or even 10 -3 M©yr _1 for 
the most intense bursts (see inset in Fig. 3). For such rates, a 10 Mj up accreting protostar mass is built up 



within ~ 10 to 10 yr. It is easily seen from eqns.(Bl) and (B3) that in that case P acc «(10-1000) xP p h, i.e. 
« 10 5 -10 dyne cm -2 , the upper values corresponding to the most intense burst events. Eqn.(B2) predicts 
that, at constant mean opacity, r oc P, i.e. a factor ~ 10-1000 increase in pressure corresponds to about 
the same increase in optical depth. As temperature increases with pressure, however, the Rosseland opacity 
increases (for P ~ 10-100) by a factor ~ 10 between ~ 2000 and ~ 4000 K (Fergusson et al. 2005). The 
optical depth at the shock front is thus likely to be of the order of r acc w 10 2 — 10 4 . Therefore, even though 
the shock energy, if the shock is supercritical, will be entirely radiated away at the shock front, this energy 
will be at least partly absorbed by the radiative precursor, since the shock occurs in a very optically-thick 
medium. At these levels, the thermal gradient is adiabatic, implying that transport by convective motions 
is more efficient than by radiative diffusion, so that the dissipated shock energy will eventually contribute 
to the proto-star/BD internal heat budget. We can try to estimate this amount of absorbed energy. 

In the optically-thick regime, the temperature at the shock front, T acc , can be crudely estimated as: 
Pace ~ Pcff Talc (e.g. Hartmann et al. 1997), which yields P acc ~ 3500-10000 K for the aforementioned 
Pacc/Pph range. The density at the shock front is thus: p acc ~ ~ 10~ 6 -10~ 5 g cm -1 . 

Using the hydrostatic equilibrium equation, VP = pW(r), where V(r) = GM/r is the gravitational 
potential at distance r from the center, and assuming AP = P* — P acc <C P*, where P* denotes the 
protostar radius and P acc the radius at the shock front, so that AP represents the depth penetration of the 
shock, one gets the following rough estimate: 

^ w MP „ 5 x io- 4 x (B5) 
P P GM 5X10 X( 106 K 10-6 } { M r J { R G > (B5j 

where AP = P acc — P p h » P acc (see above), and p ps p acc - For P acc s» 10 6 dyne cm~ 2 (see above) this 
yields AP/P ~ 0.05, while for the most intense accretion bursts, with P acc « 10 6 — 10 7 dyne cm~ 2 , one 
gets AP/P > 0.1, for ~ 10Mj up . According to these estimates, the shock, for such high accretion bursts 
can thus penetrate deep enough below the stellar surface to reach the internal adiabat, which almost reaches 
the photosphere for low-mass objects. As mentioned above, the absorbed energy should thus be efficiently 
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redistributed within the entire protostar interior, modifying its energy budget. The corresponding fraction 
of accreting energy absorbed by the protostar, compared to its binding energy E, can be estimated as 



E 



acc 




(3 - 10) T cff 

(T) 



0.03-0.10 



(B6) 



where cy is the specific heat per unit mass and (T) w w 10 7 (]^)(^) Hixf) 1 K. Larger accretion 

rates, as expected for the largest cores, will lead to even larger energy contributions. These estimates tend to 
corroborate corresponding values of a for masses issued from high burst accretion phases, as indeed suggested 
in the text for the Fu Ori events. 

As mentioned above, these estimates should be taken with great caution and detailed numerical calcu- 
lations of the accretion shock properties at the second core stage are strongly needed in order to get more 
accurate determinations. These estimates, however, seem to bring support to the results obtained in the 
calculations presented in the text, namely: 

(i) a a — > value for the lowest mass objects (mainly within the BD domain), as these objects result 
from the collapse of prestellar cores which hardly ever produce accretion rates exceeding the mean value 
cjg/G ~ 10~ 6 M©yr . This leads to accretion shocks which remain located near the photosphere level or at 
worse dissipate a negligible amount of energy within the protostar interior. 

(ii) possible values a > 0.1 for higher mass cores, thus higher mass stars, within the stellar domain, 
since the collapse of these cores can lead to accretion bursts several orders of magnitude larger than the 
aforementioned mean value. Such shocks can penetrate within a substantial fraction of the forming protostar 
and thus dissipate their energy in the inner regions, contributing non negligibly to the protostar's energy 
budget, halting or even counteracting its contraction. 

This globally suggests cold accretion conditions (a — »■ 0) for the very low mass domain and possible "hot 
accretion" events for higher masses. 
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